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Abstract 

The spectrum of eigenenergies of a quantum integrable system whose hamil- 
tonian depends on a single parameter shows degeneracies (crossings) when the 
parameter varies. We derive a semiclassical expression for the density of cross- 
ings in the plane energy-parameter, that is the number of crossings per unit of 
energy and unit of parameter, in terms of classical periodic orbits. We compare 
the results of the semiclassical formula with exact quantum calculations for two 
specific quantum integrable billiards. 
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INTRODUCTION 



The analysis of energy spectra of different physical bounded systems has always been an 
interesting subject in quantum mechanics. The energy level spacings as well as the existence 
of degeneracies (crossings) in such systems have been widely studied during recent years. The 
number of degrees of freedom, the separability of the problem and the number of free parameters 
involved, among others, are very important features which have to be taken into account when 
a given quantum spectrum is analyzed. 

Starting from Percival's ideas |I| for systems with more than one degree of freedom, two kinds 
of spectra have been distinguished. The regular spectrum, whose level spacing is characterized 
by a Poisson distribution is associated with integrable problems (the harmonic oscillator is 
an exception). The other one, the irregular spectrum closely characterized by spectral statistic 
of the Gaussian ensembles, corresponds to non-integrable systems This statistical behav- 

ior can be related to crossings and repulsions between levels with the same symmetry when 
the spectra are analyzed as a function of the parameters of the problem. In fact, integrable 
systems depending on one parameter exhibit many crossings while non-integrable systems show 
repulsions and double- hyperbola curves (avoided crossings), rather than degeneracies. So, for 
one parameter dependent systems, the presence of crossings or avoided crossings are quantum 
fingerprints of the properties of the classical dynamics. When the system depends on two param- 
eters, the surfaces of energy (that define in such cases the eigenenergies) will cross in continuum 
curves leading to diedric intersections when the system remains integrable or they will intersect 
at isolated points ("diabolical points") if the systems is not integrable (assuming that it has 
time reversal symmetry) ||. 

Distributions of avoided crossings (for one parameter dependent systems) according to prop- 
erties such that the closest approach or the mean and difference between the slopes of the 
involved levels were already established for a generic quantum system employing parametric 
random- matrix models [|7]||. The same models, extended to two and more parameters, were 
used to study distributions of diabolical points M and other relevant distributions of singular 
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On the other hand, it is well known that classical periodic orbits are essential elements 
to develop semiclassical quantization methods. Not only for integrable systems (through the 
Berry- Tabor formula Jll| ) but also for the non- integrable ones (using the Gutzwiller methods 
||12|| ) the spectral density can be formally described in terms of closed orbits of the classical 
system. Therefore, it would not be surprising that other quantum densities such as densities of 
degeneracies or densities of avoided crossings can be related to the periodic orbits. 

In the present paper, we find the density of degeneracies as an expansion in terms of the 
periodic orbits for classical integrable systems depending on one parameter. 

In section | we introduce the density of degeneracies. The semiclassical version can be written 
as a sum of a smooth part and oscillating contributions depending on the periodic orbits of the 
classical system. Section || is devoted to compute the smooth part of the density of degeneracies 
and related distributions for two specific integrable systems whose hamiltonians depend on a 
single parameter in a different functional form. We study the rectangular billiard of sides a and b 
where the parameter is the ratio fi = b/a (shape parameter) and the Aharonov-Bohm cylindrical 
billiard where the parameter is the magnetic flux, that is fi = <fi. The oscillating contributions 
are computed in Section |TTJ. Finally Section fV| is devoted to concluding remarks. We have 
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included an Appendix that contains the appropriate derivations of the semiclassical density of 
degeneracies when the quantum system has not Kramer's degeneracies (like the Aharonov-Bohm 
cylindrical billiard). 



I. THE DENSITY OF DEGENERACIES 

We consider an integrable system whose hamiltonian depends on a single parameter p. That 



is: 



H=H{I,fi) 



XI) 



— » 1 

where I = (ii, .., J n ) and h = ^ §Pid<li are the action variables. To obtain the eigenvalues, we 
can employ the E.B.K. semiclassical quantization rule 



h{ni + an/4) 



1.2) 



where n« = 0, 1, 2, ... and a, are the Maslov index [|12|]. So, we establish the quantum eigenener- 
gies E by 



E(n, fi) = H(n + a/4, p) . 



1.3) 



When we consider the eigenenergies as a function of the parameter p, the spectrum shows 
degeneracies (crossings). They occur whenever 



H(n + a/4, p) - if (n' + a/4, /i) = . 



11) 



Given n and n', this equation determines the values of the parameter p for which the eigenen- 
ergies labeled by n and n 1 are degenerated, 



/i = L(n + a/4, n' + a/4) . 
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We define the density of degeneracies p c (E,fi) as the number of crossings that occurs in the 
energy interval [E, E + dE] and in the parameter interval \p, [i + dfj] . 
Therefore, using the E.B.K. rule, we can write p c (E,fi) as follows: 



Pc (e, fi) = - 5 ( E - H ( H + «/ 4 ' - L ^ + «/ 4 ' H> + «/ 4 )) 



1.6) 



where we disregard as in Ref. || possible degeneracy factors. Employing the Poisson sumation 
formula in Eq. (|1.6|) we write, 



Pc(E,n) 



2h 



2/i 



m m' 
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\ \ d n I d n I'5 (E -H(I,i4)d(ii- L(I, ?)) exp 
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In the following, to make the expressions more handled, we assume two degrees of freedom. 
To eliminate the (^-functions, we change the integration variables as follows, 



where we have defined 



i!=E-H{I u I 2 ,n), 
£2 = p - L(h, h, I[, I 2 ) . 



Therefore 



dlidl2dl[dl' 2 



dH dL 
dh dl' 2 



dlidl 2 d^id^2 



The partial derivative dL/dI 2 is defined by the implicit equation 

H(I 1 ,I 2 ,L)-H(I[,I! 2 ,L) = 

and we obtain 



dL 
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dH 



dH\ dH 



dL 



dt 



dK ' 
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where §f-)j (ff")/') * s the partial derivative evaluated in / (I'). After integration over £1 and 
£2 follows 



fh m' 



-t — la-m + a-m 
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where 
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The first term in the Eq.( |1.13| ) (with m\ = m 2 = m[ = m 2 = 0) corresponds to replace in 
expression Eq.( |Q)| ) the quadruple sum by the quadruple integral. It is a smooth function (i.e. 
non oscillating) of E and /x, 



(Pc(E,fj)) = ^4 



2Ti Jh Ji[ 



dhdI ,m-w\ 



OJ2 U 2 
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Starting from Eq.( |1.18|) we can determine other relevant distributions. As an example, we can 
obtain the smooth part of the distribution of crossings (p c (E,fi, V)) according to the difference 
between the slopes of the levels, 



A 



V=\d.H-d.H'\ 
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To obtain such distribution (in the following DCDS), we perform the change of the integration 
variables in Eq.( 1.18Q 



(h,I[) -(Ji.V) 



Therefore 



dI x dV 
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where the tu's must be considered as a function of I±, V and /x. Leaving out the integration over 
V, follows 



(p e {E,ii,V)) = V 



dh 



h u 2 luj^d^ - uj'td^l 
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We stress that the dependence on V in Eq. (|1.22| ) is not only given by the prefactor V but also 
by the integrand and the limits of integration in the integral over I x . 

The other terms in Eq.( p..l3[) contain oscillating functions and we utilize the stationary 
phase technique to evaluate them. The conditions of stationary phase lead to the periodicity 
conditions, namely 

u>i mi 
u 2 



m 2 
m! x 
m' 2 
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Because u^s (u^'s) are positive, the sums over m fc 's (m' fc 's) in Eq.( |1.13[ ) are restricted to the 
first and the third quadrants. The numbers mi,m 2 define the topology of the periodic orbits. 
However, there could be pairs of periodic orbits that having the same topology, they are related 
to each other through time reversal transformation (non-self retracing orbits). 

For the terms mi ^ 0, m 2 ^ 0, m\ = m 2 = (mi = m 2 = 0, m\ ^ 0, m' 2 ^ 0) we perform 
first the integration over the variable I[ (Ii) and next we evaluate the second integral using the 
stationary phase approximation. Finally, we obtain 



oscl 
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h 7/2 



E 
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mi ,rri2 



mi,m 2 ^(0,0) 



m 2 



dlf 



1/2 



I[>0 



p{nj[)di[ 



where we have defined 
P(X,Y) 
S(m 1 ,m 2 ) 
9 (mi, m 2 ) 



d lx H(X) - d.H'(Y) 



uo 2 (X)uj 2 (Y) 

2n (mj° + m 2 l° 2 ) 

vr . / d 2 I 2 \ 
-sig m 2 
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dlf 



(«imi + a 2 m 2 ) 



cos (S(mi, m 2 )/h + 9 (mi, m 2 )) 

(1.25) 

(1.26) 
(1.27) 
(1.28) 
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m\ ,rri2 



2 if there are two orbits of topology mi,m 2 
1 otherwise . 
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1° = J°(mi,m2),/| = I 2 (mi,m 2 ) are the values of the actions given by Eq. (|1.23 ). On the 
other hand, for mi ^ or m 2 7^ and m'i 7^ or m' 2 7^ we obtain 



_2_ sr- x ^ A mi , m2 A m / i;m ^P (If (mi, m 2 ),I[°(m' l , m' 2 )) 

a mi,m 2 ^(0,0) mi,m^(0,0) 



/ rf2/ 2 d?I> 2 



1/2 



x cos (S(mi, m 2 )/Ti + 0(mi, m 2 )) cos (S'(m / 1 , m' 2 )/h + 9'(m' 1 , m 2 )) ■ 
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II. THE SMOOTH PART OF THE DENSITY OF CROSSINGS 



In the present section we compute the smooth part of the density of crossing starting from 
expression Eq.( |1.18| ) for two specific systems. We will consider two kinds of billiards whose 
hamiltonian depend on the parameter fi in a different way. 



A. Rectangular billiards 



At first we study the well known rectangular billiard. That is a spinless particle in a two 
dimensional rectangular box of sides a and b. The hamiltonian in terms of the action variables 
is 



7T 



(2.1) 



In this case we will consider the crossings as a function of the shape parameter /1 = b/a. We fix 
the area of the box A = ab as a constant to conserve invariant the smooth part of the density 
of states. So the hamiltonian can be written as 



I 2 

h = fill + ± 
fi 



where h 



2mAH 



. In the following we will use the hamiltonian given by Eq.(|2 



(2.2) 

and we call the 



corresponding energy as e = 2mJ ^ E . 
Taking into account that 
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<jJ 2 



fie - fi 2 l\ , 

> 

fi 



{d^h-d^til =2 
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1 l 1 l 
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and setting h = 1, we obtain, after replacing in Eq.( |1.18| ) 



fi h=V^ dl[ 
{Pc e,/i = 7 / 1 

4 Jl i=o Je - fil[ 2 



h=V^ (i 2 - r 2 )dh 

h=K Je-fil 2 



+ 



h=I[ (I[ 2 - I 2 )dh 



fill 
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that after elementary integrations leads to 



<Pc(e,AO> = ~- (2-7) 

To compare this result with the exact one obtained by quantum calculation, we integrate 
the above density of crossings over the parameter /x. In this way we establish the mean number 
of crossing in the finite interval [//i,/^] with energies between e and e + de. 

I dn \ r^ 2 

1 c (e,/^i,/^2) ) = / (p e (e, n)) dfi, 



de 




{21 



which is a linear function of e. Figure [I] shows the predictions of the preceding result for /xi — 1 
and Hi — 2 (solid line) and /X2 — 6 (dashed line) together with the exact calculations. 

Now we apply Eq.( |1.22j ) to obtain the DCDS. We define the relative difference of slopes as 



v = V/V max where V = 2 |If — If | and V max is the highest value of V for a given value of the 
energy e and the parameter /x (in the present example V max = 2e//x). By changing variables in 
Eq.( [2.6| ) — > (z — Ii\Jp,/e,v) and leaving out the integration over v, we obtain 

(Pc(e,fi,v)) = ]-g(v), (2.9) 
4 /z 



where 



g(v) = v [ Z 1 . dZ == . (2.10) 

J*=>fi z 2 ){l-{z 2 -v)){z 2 -v) 



Therefore, the joint distribution Eq.( p.9|) results factorizable as (p c (e, /-0) ^("y)- This means that 
the distribution of crossings according to the relative difference between the slopes of the levels 
is a global property of the system. It holds for any region in the plane e — fi. Equation ( f2.10| ) 
is normalized to one and it gives the fraction of crossings whose relative difference of slopes 
is in the interval [v, v + dv] and it is independent on the values of e and /x. For example, the 
distribution of the number of crossings in the intervals [0, e] and [1, /x] according to the difference 
of the slopes is 

e 2 In a , , . 
= —^9(v) • (2.11) 

That is, if we perform an histogram according to the relative difference between the slopes for 
all the crossings in the intervals [0, e] and [1, /x], the smooth part of such histogram will be given 
by Eq.( 2.10] ). Figure |2| shows the distribution Eq.( p.lOj ) together with the histogram resulting 




from the exact quantum calculation computing the relative difference between the slopes of the 
crossings that occur in the region e < 3000 and 1 < /x < 2 (~ 720000 crossings). 



B. Aharonov-Bohm cylindrical billiards 

As a second example, we consider the Aharonov-Bohm cylindrical billiard. That is a spinless 
particle confined in a two dimensional cylindrical shell of height a and radius r axially threaded 
by a confined magnetic flux <fi. The hamiltonian is (taking Ii non negative) 



7 



2mr 2 



2ma 2 



that can be rewritten as 



h = ± <pf + L 



2 ■ 



where 7 = -§-5, 



and /i 

2-7TC 



2ma 2 H 



. As in the first example we define the energy e 



(2.12) 
(2.13) 

2ma?E 



Now, we will consider the crossings as a function of the normalized magnetic flux 
set \i = (p. Taking into account that 



and setting Ti 



1, we obtain: 



1 2 = ^/e-7(/i ±//) 5 
oj 2 = 2/ 2 , 

27|A±/il , 



That is we 

(2.14) 
(2.15) 
(2.16) 



(Pc(e,/i)) =7 



/e/7 



+ 



4=o ^ - 7 /f 
Ji=4 (J( - h)dh 



in 



h=y/e/>y _ J/)^ 

/i=V 

/i=0 



2 v /e - 7 J? 



7i? 



(2.17) 



To derive Eq.( 2.17 ) it is necessary to take into account in Eq. (|1.6j ), the effect of the breaking of 
the time reversal invariance (see Appendix). After performing the integrations, Eq. (|2.17|) gives 



(Pc(e,/x)) 



2^ 



(2.18) 



This smooth density of crossings is independent on the parameter [i (the flux). Therefore the 
number of crossing between e and e + de per unit of flux is given by Eq. ( |2.18j ) . Figure [5] shows 
(p c (e, //)) given by Eq.( |2.18|) and the exact quantum calculation for a system with 7 = \. 

Starting from Eq.( p,17|) , it is not difficult to establish the DCDS. Let distinguish the crossings 
according to the relative sign between the slopes. We label by a plus sign (+) (minus sign, (— )) 
the crossings between levels with equal (different) sign of their slopes. We can discriminate in 
expression Eq.( |2.17| ) the contributions from both kind of crossings. The first and the second 
integral in the square bracket correspond to crossings with the same sign of the slopes while the 
third integral corresponds to crossings with different sign. Thus we find 



(Pc{e,fi)) 
Pc(e,/i) 

Pc(e,/i) 



p c (e,^) + ) + (p c (e,/i) ) 




(2.19) 



Although the smooth part of the density of crossings depends on the energy e and on the shape 
parameter 7, the fractions of each kind of crossings are the same for all cylinders and they are 
independent on the energies. 



To establish the DCDS we proceed as follows. We define the relative jump of the slopes as: 

-\h-l[\ , (2.20) 



for crossings with same sign of the slopes and 

v = \[{{h + I[) , (2.21) 

for crossings with different sign of the slopes. In this way v results < v < 1 for crossings 
labeled (+) while it is < v < 2 for crossings (— ). 

Thus by changing variables in Eq.( |2.17|) and leaving out the integration over v we obtain 



Mw)) = < 2 v (e) - if 2 >,>i, (2 ' 22) 



where 



, , , 1 f 1 dz 
9i{v) = t; v 



2 Jv l-( Z - V )2^T 



z 



2 



4 ^0 ^1-( W - Z )2 V ^72 

/ x 1 f 1 dz , 

9ii(v)~ = -v / , . 2.23 

4 j v -i /i _ i v _ 2 )2yr^2 



Equation ( |2.22| ) gives the mean number of crossings per unit of flux that occur in the interval 
of energy [e, e + de] and such that the relative difference between the slopes lies in the interval 
[v, v+dv] per unit of flux. For v < 1 there are two contibutions. The first integral corresponds to 
crossings with the same sign of the slopes and the second one corresponds to crossings between 
levels with different sign of the slopes. For v > 1, only crossings with different sign of the slopes 
can occur. Let us remark that the joint distribution (p c (e, /x, v)) results factorizable (as in the 
rectangular billiard) 

(p c {e,fji,v)) = (p c (e,/i)) g(v) , (2.24) 

with 

g(v) = I 9 'W + 9 'W~ I!!- 1 ' 9 (2.25) 
[ 9ii( v ) if 1 < v < 2 . v ' 

As in the rectangular billiard, distribution Eq.( |2.25| ) holds irrespective of the values of the flux 
and the energy. Therefore, the distribution of crossing below an energy e according to the value 
of v will be given by 

/ dn c (e,n,v)\ r , , 

-) =9{v) / (pc{e,fi))de 



dv / 

4 (e) 3 / 2 



o - 9(v). (2.26) 
3 ^fl 



Figure f| shows g(v) predicted by Eq.( |2.23j ) together with the exact quantum results obtained 
computing the crossings for e < 1400 between < /x < 1 (n c = 105158 crossings). We have 
discriminated the crossings according to the same sign (gj(v) + , dashed line) or different sign 
(gj + gjj, solid line) between the slopes of the levels. We remark that there are n+ = 22266 
crossings (+) and n~ = 82892 crossings (—). These results are consistent with Eq.( j2.19| ) that 
implies 

i£ = 2 - _ 

n c 2 ' 

£ = ' (2.27) 
n r 2 



III. THE OSCILLATING CONTRIBUTIONS 

In this section we calculate the oscillating contributions to the density of crossings Eq. ( 1.25| ) . 
They are determined by the periodic orbits that satisfy condition Eq.( |L23j ). 



A. Rectangular billiards 



Given the hamiltonian Eq.( |2.1| ) the periodic orbits are classified according to the topology 
(mi, m 2 ) being mi(m 2 ) the number of bounces on a side of length 6(a) before the periodic orbits 
are closed. 

Taking into account Eq.( |2J]) and the requirement Eq.( |1.23j ), we have: 



I°2 



mi 
m 2 



and 



I°(m 1 ,m 2 ) ■■ 

I°(m 1 ,m 2 ) ■■ 
Therefore from Eqs. ( |1.27| ) and ( |1.28| ) follows: 

S(mi,m 2 ) = 2n t 
6(1111, m 2 ) = 



mi 




m 2 /i 



2 

ml 



h m 2/ u e 



7T 
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(3.1) 



(3-2) 



(3.3) 
(3.4) 



On the other hand, using expression Eq.( pT3|) , we have 



in 



dl 2 l 



/x 3 e 



(3.5) 



We replace Eqs.( |3T2|),( |3.3| ) and ( |3.5| ) in Eq.( pL .25| ) and performing elemental integrations we 
obtain the oscillating contributions: 



Pc(e,n)osd = (-) 

mi ,rri2 ,4(0,0) W 

| (mj - mlii 2 ) 
2{m 2 + m 2 ,/! 2 

X COS 27T 



\ 3/4 



(m 2 +mln 2 ) 
mi 

arcsm ' 



(4mim 2 /i + nm 2 ,^ 2 



7rm; 



?(m 2 + m\n 2 ^ 



\jm\ + m 2 ,/! 2 



On the other hand, replacing Eqs. 

Pc(e> fl)osc2 = 



Q[) and 



mi,m 2 ^(0,0) mi,m 2 ^(0,0) \^/ 



in Eq. ( |1.30|) we obtain 
1 

{ml + m 2 ii 2 yl\m! 2 + m' 2 ^ 2 ) 1 ^ 



x 



/2 



(m 2 + m 2 .// 2 ) 



(m /2 + m' 2 2 /i 2 ) 



x cos 27T< 



mf „ . 7r i I 



(3.6) 



(3.7) 



Unlike the smooth part Eq.( |2.7| ), both contributions involve information about the individual 
crossings. The calculation of p c (^,p)osc2 implies a quadruple sum over m's while p c (e, p) osc i is 
a double sum. However if we consider only p c (e,fi) osc i we have a satisfactory resolution of the 
individual crossings in the plane e — /i as can be seen in Fig. [|. This figure shows a region in the 
plane e — fi where there are three crossings. In Fig. |5](a) the exact quantum levels are plotted. 
Figure ||(b) shows the oscillating part Eq.( |1.25| ) considering mi and m 2 up to 150. 



B. Aharonov-Bohm cylindrical billiards 

In this example, the periodic orbits are labeled by the number of rotations around the axis of 
the cylinder (mi) and the number of bounces on the base of the cylindrical surface (777.2) before 
the periodic orbits are closed. Taking into account that 
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mi 
m 2 
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S(m 1 , m 2 ) = 2^ 
7e 



7 



(ml + ym$) , 



d 2 h 
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We obtain (for details see the Appendix): 
Pc(e,v)osci = Y A 

mi,m 2 ^(0,0) 

^2 
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(3.9) 
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where A 
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x cos f 2nJ —(ml + 7ml) I cos (27rm]u) , 

V V7 "4/ 



2 if mi 7^ and A = 1, and 

mi 
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Pc(e,p)osc2= Y Y 
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m, 



(m 2 + 7ml) (mf + 7m 2 2 ) 



cos 



27T ( J — (mf + 7m 2 ) — J—(mf + 7m' 2 2 ) ] ] cos (27r(mi — m'^/i) 
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+ cos ( 2tt ( a —(ml + 7m 2 ) + 1/ — (m /2 + 7m' 2 2 ) — 1/4 I I cos (27r(mi + m[)p) 
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mi 
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m. 



(m 2 + 7m 2 .) (rn'i 2 + 7mf) 



cos ( 27r I J — (ml + 7m 2 ,) + * —(mf + 7m' 2 2 ) — 1/4 I I cos (27r(mi — m'-^p) 
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+ cos ( 27r ( J— (ml + 7m 2 .) — J—(m'i + 7m' 2 2 ) ] ] cos (27r(mi + m'-^p) 



7 



7 



. (3.11) 



As in the example of the previous subsection, p c (e, p) osc i Eq.( p.lOj ), gives enough information 
to determine individual crossings, as can be seen in Fig. |5|. In Fig. |^(a) is depicted a region 
of the plane e — p where the exact quantum levels show three crossings. Figure |6|(b) shows the 
oscillating contribution given by Eq.( p.lOp taking into account mi and m 2 up to 150. 

Unlike the rectangular billiard, in the present case, the dependence on the parameter p 
of the oscillating contributions is quite simple. Therefore Eqs.( p. 10 ) and ( |3.11|) can be easily 
integrated to obtain the oscillating part of the number of crossings in the energy interval [e, e+de] 
per unity of flux, 



dn r 



de 



— = / p c (e,p)dp 



(pc) dp + p c (e, p)oscidp + / p c (e, p) OSC 2dp . 



fj,=0 



(3.12) 
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As we have pointed out in section [TT], the smooth part of the density of crossings depends 
only on the energy e, so the integrated smooth part (the first integral in the right hand side of 
Eq.( |3.12| )) is given by Eq.( [2.18| ). For the oscillating contributions, the integration over p of the 
terms in Eq.( |3.10|) vanishes unless mi = (because they are proportional to sin(27rmi)/27rmi). 



Therefore we obtain: 

=1 



VT/ m2 >o m 2 




p c (e,p) oscl dp = 2[-) —m cos ( 27rm 2V^ ~ 7) • (3-13) 



In terms of the classical motion, the contribution Eq.( |3.13| ) is originated by the orbits that 



having m 1 = 0, correspond to an irrotational motion with bounces between the bases of the 
cylinder. Figure |7|(a) shows the exact quantum calculation of the integrated density of crossings 
obtained through an histogram. The sum of the contribution Eq.( |3.13| ) (up to m<i = 500) and 
the smooth part Eq.( |2.18| ) is shown in Fig. 0(b). We can see that the contribution of the 
oscillating term Eq.( |3.13|) originates the sharp peaks that are present in the exact calculation. 
By inspection of Eq.( |3.13| ), we expect a sharp peak whenever e tunes a value such that 

m 2V /i-l/8 = Z, (3.14) 

for all m 2 being / an integer. For e >> 1 this condition leads to e ~ n 2 . Therefore we conclude 
that the peaks correspond to the regions of energy where a "head-rotational band" state appears. 
Such states (the first state for a given n ) are the less affected by the change of the flux (p) 
(they have the smallest slope because they have angular momentum equal to zero) and they 
contribute to the density of crossings in a relevant way because they will cross with all the 
other states that go through (upward or downward) this region of energy. The integration of 
the terms corresponding to the other oscillating contribution Eq.( |3.11 ) also vanishes unless 



mi = m[. Taking into account these contributions, the integrated density of crossings is shown 
in fig. 0(c). 



IV. CONCLUDING REMARKS 

In the present paper we have derived a semiclassical expression for the density of degeneracies, 
p c (E, p) that occur in quantum integrable systems depending on a single parameter p. We have 
applied our results to two specific systems that depend on the parameter in a different functional 
form, obtaining a quite satisfactory agreement with the exact quantum calculation. Our results 
show that the density of crossings is strongly dependent on the derivative of the hamiltonian 
with respect to the parameter. Therefore, unlike the density of states, even the smooth part 
of the density of crossings can depend on the energy E and on the parameter p in a diverse 
way according to the functional form of H(p). In particular, the smooth part of p c for the 
rectangular billiard, where p is the ratio between its sides, results ~ E/p. Taking into account 
that the origin of the crossings between eigenenergies of integrable billiards as a function of a 
shape parameter has the same and very simple interpretation in terms of geometrical arguments 
fl~3|j , we expect that the dependence ~ E/p of the smooth part of the density of crossings, holds 
for any integrable billiard (like the elliptic billiard where p is the ratio between the axes of the 
elliptical box or the annular billiard where p is the ratio between the outer and inner radii). On 



the other hand, for the Bunimovich stadium |14| , there are numerical evidences that the smooth 
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part of the integrated density of narrowly avoided crossings (that is the number of narrowly 
avoided crossings per unit of E for a given interval of fi, being [i the shape parameter), shows a 
linear dependence on the energy E. This suggests that the integrated distribution of narrowly 
avoided crossings when /i is a shape parameter would have the same functional dependence on 
E (linear) for other irregular billiards. 

For the Aharonov-Bohm cylindrical billiard, the smooth part of the density of crossings is 
~ \I~E~ without dependence on fi = <fi (the magnetic flux). Such functional form should hold for 
other Aharonov-Bohm integrable billiards and this is the case for the Aharonov-Bohm annular 



billiard |15 l. Moreover in Ref. |15[ it is shown, by numerical computation, that the distribution 



of spacing in flux between crossings that a given level has, is Poissonian. The present work 
provides the framework to study properties of the crossing spacing distribution, employing the 
two point correlation function for the density of degeneracies. 



APPENDIX: 

The present appendix is devoted to show how to handle the sum Eq. ( |1.6| ) for the Aharonov- 
Bohm cylindrical billiard. This system having one degree of freedom of rotation, it has not 
Kramer's degeneracies. That is, conjugate states by time reversal transformation are non de- 
generated. 

At first, let us consider the billiard without flux. The eigenenergies are 

%i,n 2 ) =7^1 + (n 2 + l) 2 , (Al) 

where n\ is the quantum number related to the rotation around the axis of the cylinder 
(ni = 0, ±1,±2, ...) and n 2 is the quantum number related to libration motion parallel to 
axis (n 2 = 0,1,2,3). The system has Kramers's degeneracy. A quantum eigenstate (ni,n 2 ) 
and its conjugate (— ni,n 2 ) have the same eigenenergy if n\ ^ 0. To take into account these 
degeneracies in the sum Eq.( |1.6| ) over ni, n! x > we must include a factor e(ni, n'i) in each term 
such that 

f 4 if m ^ and n[ ^ 
e(m, ni) = < 2 if (m ^ and n[ = 0) or (m ^ and n[ ^ 0) (A2) 
[ if ni = and n\ = . 

When the cylinder is threaded by a magnetic flux 0, the eigenenergies of a conjugate pairs 
of eigenstates are 

h{n u n 2 , 0) = 7(ni - 0) 2 + (n 2 + l) 2 , 
h(-n 1 ,n 2 , <p) = l{-n x - <p) 2 + (n 2 + l) 2 , 

= 7 (m + 0) 2 + (n 2 + l) 2 . (A3) 

Therefore, conjugated pairs of eigenstates are no longer degenerated. In such a situation, to 
preserve the sum Eq.( L6) over ni,n'i > 0, we must distinguish two kinds of states according 



to the dependence on the flux that their eigenenergies have. Moreover, Eq.( |1.5| ) which deter- 
mines the values of the parameter for which the crossings occur, changes according to the same 
dependence. Thus, Eq.( |1.6| ) for the cylindrical billiard means 
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Pc( e > <P) = \ 1212 S(e - h(ni + 0, n 2 ))5(<p - L x {n u n 2 , n[, n' 2 )) 

n n' 

+ J E 12 S(e-h(Tn -0,n 2 ))5(0-L 2 (ni,n 2 ,n / 1 ,n / 2 )) 

n^(0,n 2 ) nY(0,ng 

+ 12 S i e - h ( n i -0>^2))5(0-L 3 (n 1 ,n2,n' 1 ,n 2 )) 

+ ^ Z Z 5 ( e ~ /l ( n i + ^' n 2))5(0-^3(^i,n 2 ,n , 1 ,n 2 )) , (A4) 

n^(0,n2) ft' 

where Li, L 2 and L 3 are the values of the flux determined by the roots of 

h(m + 0, n 2 ) - h(n[ + 0, ri 2 ) = , 

/i(ni — 0, n 2 ) — /i(ni — 0, n' 2 ) = , (A5) 
h{nx - 0, n 2 ) - ^(^i + 0, n 2 ) = , 

respectively. 

Now we employ the Poisson's formula. The first sum in Eq.( |A4|) that takes into account the 
crossings between energy levels that increase as a function of (this term includes the crossings 
between levels {n\ = 0,n 2 ) and (n[ = 0,n 2 )) can be written as, 

\ 12 12 L L 5(e-/i(/ 1 + 0,J 2 ))5(0-L(/ 1 + 0,/ 2 ,J{ + 0,^) 

2 mfX 2m ^ - ( °' 0) 7/ '- (0 ' 0) 

x exp (i2n(m ■ I + rh' ■ /')) , (A6) 

where we have taken into account that ot\ = and a 2 = 4. To elimate the 5 functions we employ 
Eqs. (1.7-1. 18) together with the first equation of Eq.( |A5| ). We obtain 



m m' — i— 



exp (%2-K(mxh + m 2 I 2 + w! x I' x + m 2 I 2 )) • (A7) 



2^ e - 7 (/i + 0)V e -7(/{ + 
Now, we redefine Ji {![) = l\ + (/{ + 0). After replacing it follow 

; exp (z27r(mi(/i - 0) + m 2 J 2 + m'^/J - 0) + m 2 I 2 )) , (A8) 



1/1-/(1 



where we have omitted the bar in the integration variables. In the same way, the second sum 
in Eq.( ~KA ) (using the second equation of Eq.( |A5| )) and the third and fourth (using the third 



equation of Eq.( [Ao]) ) are 

mm' — ^ i 



|J — l'\ 

x — 1 ? y_ exp (i27r(mi(/i + 0) + m 2 / 2 + mM + 0) + m 2 / 2 )) , (A9) 
2 Je - jlfJe - Tif 



1K 



7 

m m 

(h + A 



2 z ^^Jh> ( f>Ji' 1 >- ( t> 



2i 



exp (i2n(nii(Ii - <fi) + m 2 h + m^lj + 4>) + m^)) > (A10) 



+) = ?££/, , /, .dfidil 



2 ; ~ Jh>-6 Jl'>q 
mm' — ^ i — v 

(A + 4) 



2/ 



exp (^(mi^ + 0) + m 2 / 2 + mi (/J - <p) + m 2 J 2 )) . (All) 



Equation ( |A9|) takes into account the crossings between energy levels that decrease as a function 
of 4> while Eqs.( |A10|) and ( All|) correspond to crossings between levels that increase and 



decrease as a function of 0. For the smooth part, we set mi = m 2 = m 1 = m 2 = in the sum 
of Eqs. ( |A8| - |A11| ) and taking into account that the lower limit in the integration over li can be 



taken as (because e >> 1 implies <fi — > + and —<fi — > _ ) we obtain Eq.( |2.17| ). 

To obtain the first contribution to the oscillating part p c (e,0) O sci Eq.( |1.25|) , we set mi = 
m 2 = and we integrate over Ii. Then, we approximate the integration over I[ using the 
stationary phase thechnique obtaining Eq.( |3.10|) . For the second contribution p c (e,0)osc2 we 



evalute both integrals (over I\ and 7 2 ) in the phase stationary condition. The resulting Eq.( ^.11 ) 



is long but straightforward. However we will remark some steps. After replacing the phase 
stationary condition for the actions and reducing the sums over and to the positive 
quadrants we will have in Eqs.( |A8| ) and ( |A9| ) terms like (we write the prefactors as A ), 

2 A cos (S - 2vrmi0) - vr/4 cos (S' - 27rm' 1 (() - 7r/4) , 

2A cos (S + 27rmi0) - vr/4 cos (S' + 2nm' 1 4> - vr/4) , (A12) 

respectively. These terms can be appropriately combined using trigonometric identities in the 
form 

A [cos (S — S') cos (27r0(mi — m'i)) + cos (S + S' + 7r/2) cos (27r0(mi + m^))] . (A13) 
In the same way the terms of Eqs.( |A10|) and ( |A11|) result 

2B cos (S - 2vrmi0 - vr/4) cos (S' + 2irm' 1 <f) - vr/4) , 

2B cos (S + 2vrmi0 - tt /4) cos (S' - 2nm' 1 4> - vr/4) . ( A14) 

Therefore the sum of these terms is 

B [cos (S — 5"') cos (27r0(mi + m^)) + cos (S + S' + 7r/2) cos (27r0(mi — m^))] , (A15) 



and Eq.( p.ll|) follows. Let us remark that Eqs.( [A12j ) and ( |A14|) are a consequence of the 



breaking of the time reversal symmetry when the magnetic flux is present. In such a case, twin 
classical orbits (related to each othe by the time reversal transformation) that have the same 
action when the systems has the time reversal symmetry split their actions when the flux is 
applied Jl6| . 
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FIG. 1. Integrated density of crossings dn c /de as a function of the energy e for the rectangular 
billiard. The straight lines correspond to the smooth part (dn c /de) = | hi(^) for ^2 = 6 (dashed line) 
and ^2 = 2 (solid line). In both cases /ii = 1. 
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FIG. 2. Distribution of crossings g(v) according to the relative difference between the slopes v for 
the rectangular billiard. The result predicted for the smooth part of density of crossings corresponds 
to the dashed line. 
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FIG. 3. Integrated density of crossings dn c /de as a function of energy e for the Aharonov-Bohm 
cylindrical billiard with 7 = 4y. The solid smooth line corresponds to the smooth part (dn c /de) = 2. 
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FIG. 4. Distribution of crossings g(v) according to the relative difference between the slopes v. a) 
Exact quantum results, b) Results predicted for the smooth part of the density of crossings. The dashed 
line curve corresponds to crossings between levels with equal sign of their slopes, while the solid line 
corresponds to crossings between levels with different sign. 
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(a) 




FIG. 5. a) Detail of three exact quantum eigenergies for the rectangular billiard as a function of // 
where three crossings can be observed, b) Density plot of the oscillating part p Coscl (see III A in the 
text) of the density of crossings in the same region of the plane e — [i . 
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(a) 





FIG. 6. a) Detail of three exact quantum eigenergies for the Aharonov- Bohm cylindrical billiard 
as a function of the flux cj) where three crossings can be observed, b) Density plot of the oscillating 
part p Coacl (see IIIBj in the text) of the density of crossings in the same region of the plane e — <fi. 
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FIG. 7. a)Integrated density of crossings dn c /de as a function of energy e for the Aharonov-Bohm 
cylindrical billiard with 7 = b) The same density predicted by the sum of the smooth part (p c ) 
and the contribution originated by p Coscl ■ c) The same density predicted by (p c ) + p Coscl + Pc oac2 ( see 



IHB in the text). 
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